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$_i ■ values of BR(b — > s#) up to ~ 10% or more are compatible with the constraints coming from 
— > the CLEO experimental results on BR(B — > X s 7) and BR(B — ► X s (j>). Such large values can 
reconcile the predictions of both the semileptonic branching ratio and the charm counting with 
the present experimental results. We also discuss a supersymmetric model with gluino-mediated 
flavour violations, which can account for such large values of BR(b — > sg). 



ABSTRACT 

We consider the possibility of having a large branching ratio for the decay b —>■ sg coming 
from an enhanced Wilson coefficient of the chromomagnetic dipole operator. We show that 
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1 Introduction 



The world average of all measurements of the 5-meson semileptonic branching ratio is p| 



BRfl = BR(B -> Xcv e ) = (10.4 ± 0.4)% 



o . 



This result is considerably smaller than the theoretical prediction in the parton model, where 
BRsl ~ 13-15% [[|. Furthermore non-perturbative corrections cannot reduce the pre- 

dicted BRsl below 12.5% ||. However it has recently been found that charm mass corrections 
to T(6 — > ccs) are large and can further reduce the theoretical prediction for BRsl 0: 

j (12.0 ± 0.7 ± 0.5±?"i ± 0.2)% on-shell scheme . . 

BRsl - | (n 3 ± Q Q ± J +o.9 ± Q 2 ^ % m gcheme {*) 

Here the errors correspond to the uncertainties on mj, a s (Mz), the renormalization scale \x and 
the other parameters present in the calculation. It may now seem that the discrepancy in the B 
semileptonic branching ratio has essentially disappeared. However, as we lower the theoretical 
prediction for BRsl by increasing T(b —>■ ccs), we simultaneously increase the prediction for n c , 
the average number of charm hadrons produced per B decay. Indeed the theoretical prediction 
for n c obtained from eq. (|2|) is 

f 1.24 ± 0.05 ± 0.01 on-shell scheme , , 

Uc ~ \ 1.30 ± 0.03 ± 0.03 ± 0.01 MS scheme ^ 

where the first error corresponds to the quark mass uncertainties, the second to the a s variation 
in MS and the last to all the other uncertainties. This prediction is larger than the present 
world average of the n c measurements [JTJ] 

n e c xp = 1.117 ±0.046 . (4) 

The real problem is therefore to understand the presence of a low BRsl together with a low 
value of n c . The theoretical difficulty is summarized in fig. [l|. 

It could well be that this is just an experimental problem; measurements of n c rely on 
model assumptions about charm hadron production which have been questioned in ref. ||. 
Recently it has been proposed || that "spectator effects" can be at the origin of the problem. 
Although these effects cannot be fully computed without a knowledge of the hadronic param- 
eters, potentially they can bring BR SL close to the experimental value. They also result in a 
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simultaneous, although modest, decrease of n c 0. On the other hand, if "spectator effects" are 
responsible for the observed low value of the ratio t{K^)/t{B^) = 0.76 ± 0.05, they then tend 
to increase the prediction for BRsl, making the problem even more pressing. Recently it has 
also been suggested J7| that the data indicate the presence of 1/ tuq corrections in non-leptonic 
decays, which are not present in HQET. This has been attributed to a possible failure of the 
operator product expansion near the physical cut and a corresponding violation of the local 
quark-hadron duality. In particular, using the phenomenological recipe of replacing the quark 
mass by the decaying hadron mass in the m 5 factor in front of all non-leptonic widths, BRsl, 
n c , and r{Kb) /r{Bd) can be reconciled with their measured values 0. 

In this paper we discuss the effect of new physics on BRsl and n c . The possibility that large 
contributions to T(b — > sg) can eliminate the disagreement between experiments and Standard 
Model predictions was first proposed in ref. ||, and then studied in much further detail in ref. 
||. Our goal here is to refine previous analyses and consider new constraints. 

In sect, ^we perform a model-independent analysis of T(b — *> sg) and the possible constraints 
coming from T(B — > X s 7) and T(B — > X s <p), extending the results of ref. [f|. In particular we 
show that T(B — > X s <p), in spite of its sensitivity to penguin operators, provides a poor probe of 
New Physics effects. In sect. |3|we discuss a supersymmetric model with gluino-mediated flavour 
violations, first suggested by Kagan ||, which can explain a large enhancement of T(b — > sg). 

2 Model-Independent Analysis 

The effective Hamiltonian relevant for AB = 1 decays is given by [ |T0[ | 

r AB=l _ Gf 



TT' 



v; b v cs (dMQSM + cmq c M) + (5) 
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V: b V us (Ci(ai)Q?(aO + CMQ U M) - VZV u y EC i Qi)Q i Qi) 

i=3 

where V^- are the CKM matrix elements and Cj(/x) are the Wilson coefficients evaluated at the 
scale n of order m&. The dimension-six local operator basis Qi is given by 

Qi = (s a qf3)v-A(qpb a ) V - A 
Q\ = {s a q a )v-A{(ipbf})v-A 
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Parameter 


Value 


\v*v tb \ 2 /\v cb \ 2 

\ ZS lu 1 / 1 1 


0.95 


a s (M z ) 


0.117 


m t (GeV) 


174 


m b (GeV) 


5 


m s (MeV) 


500 


m c /m b 


0.316 


(GeV) 


1.019 


9'i (GeV 4 ) 


0.0586 



Table 1: Central values of the input parameters used in the analysis 



Q3,5 — (Sgba)v-A ^(gggg) Vt^ 

1 

Qifi = {Sgbp)y-A y^fega VtA 

3 

^7,9 = -{s a b a )v-A X! e q(.Wp)v±A 
1 q 

3 

^8,10 = ^(SqMv-A e qi.Wa)v±A 
1 q 

r\ Us — uv ,a l /^yA 

- Y^ m & S a 0r V+A t a/3°/3 G > 

g 

Here (gi^) (17394) denotes current-current products, V ± A indicates the chiral structure and 
a, /? are colour indices. Moreover, p s (e) is the strong (electromagnetic) coupling constant, G A U 
(F^) is the gluon (photon) field strength, and the t A are the SU (N) colour matrices normalized 
so that Tr{t A t B ) = 5 AB /2. 

This Hamiltonian is known at the next-to-leading order (NLO), as far as one separately 
considers the current-current, gluon and photon penguin operators Q1-Q10 on the one hand 
llj and the magnetic-type operators Qn-Qu on the other hand Hl2] . Unfortunately the mixing 



between these two classes of operators is at present known only at the leading order (LO) JT3J . 
The branching ratio of the decay b — > sg is given by |T3| 

BR(b - sg) = BR ^Yf J ^(/x)| 2 , (7) 

\V cb \ 2 irg{m c /m b ) 

where BR e $l is the measured semileptonic branching ratio, the phase-space factor g(z) is given 
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in the Appendix, \i = 0{m^) and CyJ is the renormalization-scheme invariant coefficient^ 
introduced in ref. [15]. With the input values given in table [l|, the LO Standard Model 
prediction is 

BR SM (b -> sg) = (2.3 ±0.6) x 10~ 3 , (8) 

where the error is mainly due to the variation of /x between and and to the uncertainty 
on a a (Mz). Inclusion of the known part of the NLO corrections reduces the central value and 
the /j dependence, but introduces a significant scheme dependence [14]. 

In the effective Hamiltonian approach, physics beyond the weak scale affects only the initial 
conditions of the Wilson coefficients, namely Ci(Mw)- Therefore these coefficients can be used 
to parametrize new physics effects without referring to a specific model. In the case at hand, we 
assume that the initial condition of the Wilson coefficient Cn(M w ) is an independent variable 
and define 

r„ 



Cii(M, 



w i 



w 



where 



9i{m 2 t/ m w) > 
x(—x 2 + 5x + 2) 

4(x- l) 3 



(9) 



(10) 



3x 2 log x 
2(x- l) 4 



In terms of r„, the branching ratio of b — >• sg is given by 



BR(b -> sg) = 2.73 x 10~ 2 (0.15 + 0.14 r g 



'.ID 



assuming central values for all the relevant parameters. 

An enhanced value of BR(b — > sg) affects the relation between n c and BRsl- We find 

2 - (2 + R T + R ud + 2R 4 )BR SL - 2BR(b -> sg) , (12) 



n. 



where 16 



R T 
Rud 
R4 



F(b — > ctu) 
Yib 

r(6 



0.25 



cev) 
cud 1 ] 



4.0 ±0.4 
1 (0 — > cez/j 

r(fe — > no charm) 



r(6 — > cev) 



0.25 ±0.10. 



(13) 



Using calculations in the 't Hooft-Vcltman scheme, as we do, Cl{{ 2 = C(i v 12 . Therefore we denote them 
Cn,i2 in the following. 
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BR 

SL 

Figure 1: Correlation between the semileptonic branching ratio BRsl and the number of charms 
per B decay n c . The solid band is the Standard Model prediction, including the theoretical 
uncertainties, while the dashed one is obtained assuming a BR{b — > sg) = 9%. The experimental 
data point is also shown. 



Here T(b — > no charm) is the sum of all B decay widths into charmless final states different from 
sg. Notice that in eq. (|12|) we have eliminated the dependence on T(b — > ccs), which is the main 
source of theoretical uncertainty. If we require that the experimental central values of BRsl and 
n c , eqs. (P and (§), lie within the theoretical uncertainty, we need 7% < BR(b — > sg) < 11%. 
This corresponds to 10.2 < r g < 13.1 or —12.3 < r g < —15.2. For instance the effect of 
BR(b — > sg) = 9% in eq. (|T2|) is shown in fig. [1] (see the dashed band). 

We now want to perform a model-independent analysis of the constraints on r g . The first 
constraint comes from the observation of the inclusive decay B — > Xg'-f. The LO expression for 
BR(B 



X sl ) is|||| 
BR(B 



Xsl) 



BR71 



V*V tb \ 2 3a e 



\cW{»)\ 2 . 



(14) 



| V cb | 2 2ng(m c /m b ) 

It is well known that this branching ratio has very large QCD corrections. It also has a 
significant theoretical uncertainty, mainly due to the scale dependence in the Wilson coefficient, 
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to be taken into account. Also in this case, the known NLO terms decrease the central value 
and the /x dependence, but introduce a sizeable scheme dependence, to be cancelled by the 
unknown terms. 

The dependence on r g in eq. (|T4"| ) comes from the operator mixing in the QCD renormal- 
ization group equations. Analogously to eq. ([5J), we define 

C 12 (M W ) 



C? 2 M {M W ) 



where 



c s 12 m {m w ) = AKK), 



(15) 



(16) 



fi(x) 



X ( 



-8x 2 -5x + 7) x 2 (3a; - 2) log x 



i2(x- iy 



2(x- l) 4 



We can now rewrite eq. (14) in terms of r g and Try CIS 

BR(B -> X s7 ) = 7.10 x 10~ 4 (0.017r 9 + 0.273 r 7 + 0.313) 2 , 



(17) 



where the theoretical error is not shown. This formula gives a model-independent parametriza- 
tion of BR(B — > -XVy). Figure |2| shows the ranges of r g and r 7 allowed by the measurement 
BR{B — > X s 7) = (2.32 ±0.67) x 10~ 4 |T7|]. The existence of two different bands corresponds to 
the overall sign ambiguity of G^ 2 (/•*) in eq. (|HP- The size of the bands accounts for both the 
experimental error at the la level and a (linearly added) theoretical error of 30% JT5fl . Figure 
shows that very large enhancements of Cn(Mw) are allowed, since they can be compensated 
by appropriate, and relatively small, changes of Cui^Mw)- A large values of Cn(Mw) can 
result in a modification of the photon energy spectrum in B — > and could be tested in the 
future if the threshold in the observed E 1 is appropriately lowered^ The solution of the 
BRsL-n c problem in terms of T(b — > sg) singles out the four disconnected regions in the r g -r^ 
plane shown in fig. |2|. 



The extraction of analogous constraints from the measured upper limit [19 



BR(B -> X s (f>) < (1.1 - 2.2) x 10" 



(90% CL) 



is more involved. Available theoretical estimates |2D], |2T| rely on some assumptions which are 
questionable. In particular, hadronic matrix elements are calculated by combining perturbation 
2 We thank A. Ali for pointing out this possible experimental test. 
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-4 -2 2 4 

Figure 2: Correlation between the ratios r g andr 1; defined in eqs. (Qj and The scale on the 

right shows the corresponding values of BR{b — > sg). The regions outside the two oblique bands 
are excluded by the measured BR(B — > X s j) . The region above the dashed line is excluded by 
the upper limit on the BR(B — > X s <p), assuming £ = 1/3. However, the theoretical error is large 
and not completely under control in this case, so that the limit is plotted only for illustrative 
purpose. Finally the two oblated regions delimit the ranges allowed in the supersymmetric 
model with LL (dotted line) and LR (solid line) flavour mixing obtained for m g > 200 GeV and 
rh > 100 GeV. 

theory with the factorization method. Moreover they made a strong assumption on the quark 
momenta distribution inside the <fi meson (p s = p s = p^/2). These hypotheses result in a 
pure two-body decay b — > s<p, while for example string fragmentation models indicate high 
multiplicity final states for the b —>■ sg decay p2| . The same assumptions are also used to 
extract the experimental limit from the measurements. Thus, constraints coming from this 
limit have large theoretical uncertainties and should be taken with caution. 

Nonetheless, in the following, we calculate the constraints in the r g -r 1 plane, coming from 
the limit on BRiB — > X s (p), under the same assumptions made in refs. |Z0|, [HJ], which allow a 
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straightforward calculation of the branching ratio. In terms of r g -r 1} the following expression 
is obtained: 



BR(B —> X s (f)) = a^C^mb)] 2 + a 2 a s (m b )Re(C^(m b ))r g + a 3 a e Re(C^(m b ))r. 

+ a & a s (m b )a e r g r 1 , 



+a A a s (m b ) 2 r 2 + a 5 a 2 e r 2 + a 6 a s (m b )a e r g r 7 , (19) 



where the coefficients are functions of the masses m b , M^, m s . They are explicitly given in 
the Appendix. 

Beside C\\ and C l2 , eq. (|I9|) depends on the coefficient C^. It appears when the matrix 
elements of the relevant penguin and electropenguin operators are evaluated using the factor- 
ization method and it is given by [|0J 

CM = ct ff ^) + ci^^) + ct ff ^) + i{ct ff ^) + ci ff ^) + ci ff ^)) 

~\ W f ^) + C 9 e// (/,) + CSf'O*) + £ (Cl fS ^) + C e 9 ff (v) + Cg'M)! , (20) 



where £ = 1/N. Here iV is the number of colours, taken as a free parameter to account for 
possible deviations from the factorization result. Since the magnetic operators contribute to 
eq. ( TT9| ) through 0(a s ) (0(a e )) matrix elements, it is mandatory to calculate at NLO. 
The effective coefficients^ C| ^ — Cffl are defined so that they include the NLO contributions 



coming from the matrix elements at 0(a s ) (0(a e )), calculated in perturbation theory [20, 23 



With such a definition, Cf — C{q become complex functions. We note that the calculations 



of these coefficients in refs. [^(J, [23[] are not complete. A full NLO perturbative calculation 
of the matrix elements would allow for a complete cancellation (at the considered order) of 
the scheme dependence of the NLO Wilson coefficients that appear in C^. In our case, we use 
NLO coefficients C\ — C\q calculated in the 't Hooft-Veltman scheme, including only the largest 



contribution to the NLO matrix elements coming from Q 2 , as in ref. [20]. We have estimated 
the residual scheme dependence (~ 15%), by comparing our effective coefficients to those given 
in ref. [20]. We consider it negligible, compared with the other uncertainties present in the 
calculation. 

The presence of in eq. (|19]) could make our analysis much more involved, since this 

coefficient depends on several initial conditions other than Cn(Mw) and Ci2(Mw)- However 

3 Notice that the denomination "effective coefficient" refers to completely different definitions in the case of 
C 3 - Cio and C\\ - Ci 2 . 
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Figure 3: The coefficient \C$\ as a function of C^{Mw)-C\q{Mw), normalized to the Standard 
Model. The eight initial conditions are parametrized in terms of the four contributions coming 
from the box and the g, Z , 7 penguin diagrams. 

C<f> is quite insensitive to the initial conditions of the penguin and electropenguin operators 
Q3 ~ Qio- 111 fig- H we show |C^| as a function of these initial conditions, parametrized in 
terms of the contributions of the gluon, Z, 7 penguin and the box diagrams, normalized to 
the Standard Model. The reason of this weak dependence is twofold. On the one hand, the 
largest terms in come from the gluon penguin operators Q3 — Qq, which are insensitive to 
their initial conditions because the large mixing with Q\ — Q 2 dominates their renormalization 
group evolution down to /x = 0{m b ). On the other hand, the effective coefficients C,- get a 
large NLO contribution from the 0(a s ) matrix elements, which do not depend on the initial 
conditions of the Wilson coefficients. This implies that BR(B — > X s cf)) in general gives a poor 
probe of new physics effects, contrary to what is often stated in the literature. 

Therefore, it is reasonable to assume as a constant in the analysis, so that the BR(B — > 
X s <f)) could be used to put constraints in the r g — r 7 plane. Unfortunately this constraint 
suffers from large theoretical uncertainties. For example the predicted branching ratio changes 
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by a factor of 2 as we vary £ between 1/2 and 1/3, which is a popular way to account for the 
uncertainty of the factorization method |2(| . Moreover the assumption on the decay kinematics 
is not under control and the related uncertainty is not quantifiable. Just for illustrative purposes 
we show in fig. § the constraint coming from eq. fljjD , assuming £ = 1/3. For larger values 
of £, the bound on r g from B — > X s (j) becomes more stringent, disfavouring even further the 
solutions of the BRsL~n c problems corresponding to positive r g . However if we assign an overall 
uncertainty of a factor of 3 in the prediction of the branching ratio, any significant bound on 
r g disappears. 



An enhanced b — > sg decay rate also affects the exclusive decays B — > Kn ||21|| . We estimate 
BR(B~ — > K°7T~) ~ 10~ 5 (1 + 0.1 r g ) 2 , which corresponds to an effect of order 1 for \r g \ ~ 10. 
Again the theoretical uncertainty of this prediction is large and not under control, so we prefer 
not to show this constraint in fig. 0. 

We have seen how present constraints allow large enhancements of BR(b — ► sg) . Measure- 
ments of BR(B — > Xg'-f) require, however, a precise correlation between r g and r 7 . The main 
difficulty to solve the BRsL^n c problem in terms of new physics is to explain this correlation. 
We now turn to a discussion of models in which this is possible. 



3 An Illustrative Model 



The suggestion that anomalously large BR(b — > sg) can explain a reduction of BR$l was first 
made in ref. ||, where it was assumed that the new effective bsg interaction is mediated by 
virtual charged Higgs boson exchange. This possibility is now ruled out by a combination of 
the constraints from B — > Xrv, B — > X^, and Z — > bb. The main difficulty of models where 
the bsg interaction arises from charged-Higgs exchange, shared by most other models with 
weakly- interacting new particles, is that generically r g ~ r 7 . Constraints from BR(B — > X s ~/) 
then allow only small enhancements of r g . 

The possibility that flavour-changing quark-squark-gluino interactions can generate large 
coefficients for the chromomagnetic operator was first suggested in ref. 0. Now the loop 
generating On has a large Casimir factor, which is not present in the loop generating Ou, and 
one obtains r g /r y ~ 5-7. This allows a significant enhancement of b —>■ sg, especially if we 
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consider solutions with negative r 7 . 

We have in mind the case in which the squark mass matrix is not diagonal in the quark 
mass eigenbasis. This situation is generic in supersymmetric models derived from supergravity 
with non-minimal Kahler metric. For simplicity we consider separately two possibilities. In the 
first case, we assume that flavour non-diagonal entries of the squark mass matrix appear only 
in the left sector and there is no left-right squark mixing. The new contributions to the Wilson 
coefficients Cn and C12, evaluated at the scale of supersymmetric particle masses, are 



Cn = ^j:u; b U ls g 2 (ml/rh*) (21) 
m g i=l 

. s x(-llx 2 + 40x + 19) x 2 (x - 9) logx , x 

= 36(x-l)3 + 6(x-\r (22) 

C12 = —2 E U? b U ls f 2 (m 2 /m 2 ) (23) 

2x(— 2x 2 — 5x + 1) 4x 3 logx 

h {X) = 27(x-l)3 + 9(^1? (24) 



v/2' 



V£V U G F 



(25) 



Here is the gluino mass and in 2 are the eigenvalues of the down-squark squared-mass matrix, 
which is diagonalized by the 3x3 unitary matrix U. If rhi are nearly degenerate, eqs. ([H]) and 
can be Taylor-expanded around the common squark mass m: 

Cn = Z^f±g 3 (m 2 /m 2 ) (26) 

. , (x 2 + 172x + 19) x(x 2 - 15a; - 18) logx . . 

^ 3(X) = 36(x - 1)4 + 6(^1? (27) 

^12 = Z 6 -^Hm 2 /m 2 ) (28) 

2(-17x 2 -8x + l) 4a: 2 (x + 3)logx 
/3(X) " 27(x-l)4 + 9(x - l) 5 (29) 

Here 5m 2 LSL = J2 U* b Ui S fh 2 corresponds to the flavour non-diagonal mass insertion. 

In the second case we consider, the down-squark mass matrices for the left and right sectors 
are both diagonal, but there are flavour non-diagonal left-right (LR) mixing terms. Since these 
terms are not generated by supersymmetry breaking, in the limit of vanishing Yukawa coupling 
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constants, we will assume here that they are proportional to the corresponding Yukawa coupling. 
In this case, the mass- insert ion approximation is always adequate and the new contributions 
to the Wilson coefficients C\\ and C\2 are 



C 



m^m 6 y 
2(x + 11) (-x 2 + 16a; + 9) logx 



C 



3(x-l) 3 3(x-l) z 



= ^^%/ 4 K/m 2 ) 
4(5x + 1) 8x(x + 2) logx 



12 



(30) 

(31) 
(32) 
(33) 



9(x-l) 3 9(a;-l) 4 
The two cases here considered are the simplest because they generate only operators with the 
same chiral structure as in the Standard Model. They also describe the generic features of more 
complicated squark mass matrices. 

By varying the relevant supersymmetric parameters under the requirement that all squark 
masses are larger than 100 GeV and m g > 200 GeV, we find that the new interactions can 
generate values of r g and r 7 within the regions illustrated in fig. ^|. The prediction has a 
strong correlation in the r 9 -r 7 plane because the ratio r g /r^ depends on the supersymmetric 
parameters only weakly, through the loop functions. Thus, in spite of its generic features, the 
model can make a rather precise prediction of this ratio. It is interesting that it is indeed 
possible to reach a region in which b —>■ sg has the correct enhancement to solve the BRsi^n c 
puzzle. The solution always corresponds to the case in which the signs of both C\\ and C\% 
are opposite to the Standard Model results. In the case of left-left (LL) mixings this region is 
achieved when the mixing between the s and b squark is maximal^ and when the lightest squark 
massf] is around 100 GeV and m g around 200 GeV. This means that supersymmetry could be 
soon discovered at the Tevatron, but no light squark has to be expected at LEP2. Similar 
conclusions apply to the case of LR mixing. However if the LR mixing were not proportional 
to m&, the same effects could be obtained for much heavier squarks and gluinos, although 
the existence of colour-breaking minima could impose significant constraints on the parameter 

4 Indeed the mixing is so large that the mass-insertion approximation is not valid. The region in fig. || was 
obtained by using the complete expressions in eqs. (21) and (p3|). 

5 Here we refer to the mass of the lightest squark, which can be considerably lighter than the other squarks, 
because of the large mixing. Tevatron bounds on squark masses |^5| do not directly apply here, since they are 
obtained under the assumption of three families of degenerate squarks. 
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space. Finally notice that, in the case of a very light gluino (jrig ~ 1 GeV) the ratio r 9 /r 7 
can become much larger than that shown in fig. ^ especially for LR mixings. In this case it 
is possible to have r g as small as —10 for positive values of r 7 , and almost reach the shaded 
region in fig. 2 corresponding to positive r 7 and negative r g . 

Given a definite model of flavour violation, we can consider more tests on its consistency than 
in the case of the model-independent analysis of sec. 0. We now assume that gluino-mediated 
flavour violations occur only between the second and third generation of the down quark-squark 
sector and consider different processes with \AB\ ^ 0, which can be of experimental interest. 

The \AB\ = 2 transitions will affect the B s -B s mixing. For LL mass insertions, the new 
contribution to Am Ba is []2l| 



4„ 2 (i; 



Am Ba = -B B J B m Bs ^ 



Sm 



m 2 



G(m 2 Jm 2 ) 



(34) 



G{x) 



216 



(-2x 3 + 27a; 2 + 144x + 11) 2x(51x + 39) logx 



(35) 



(x-1) 4 (x-1) 5 
This contribution is much larger than the Standard Model one, making the discovery of B s -B s 
mixing even more arduous. However, for x ~ 2.4, the function G(x) has a zero and therefore, 
for particular values of the squark and gluino masses, we can reduce the total contribution 
to Am Bs . Nevertheless, the generic prediction of the model is that B s -B s mixing will not be 
discovered soon. 

The | AS I = 1 transitions can induce the flavour- changing decay Z — > bs. For LL transitions 



BR{Z -> bs) 



aa mz 



18tt 2 sin 2 6 W cos 2 9 



w ' 

l-x 



7 \\ 



2 
3 



1 - - sin 2 9 



U; b U ls F(m 2 ,mlm 2 z ) 



r m 2 (x + y) + m?(l — x — y) — m 2 z xy 



dx I dy log 
o Jo 



(36) 



(37) 



m 2 (l — x) + m~x 

Assuming the validity of the mass insertion approximation, eq. (BO) can be simply expressed 



in terms of r g . For m 2 = m| > m 2 z , we find 



BR(Z ^bs) = 5 x 1(T 8 r 2 



(38) 



Unfortunately, observation of this decay mode represents a real experimental challenge. LEP1 
can at best reach a sensitivity on BRiZ — > bs) of about 10~ 3 |27] . 
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Finally we consider FCNC decay processes of B mesons, such as B — > X B vv and B 



X s £ + £~ . It is known [gj] that gluino- mediated interactions do not affect the Z penguin operator. 
Indeed, the effective b-s-Z vertex turns out to be proportional to the 6-S-7 vertex. QED 
gauge invariance then implies that the Z penguin is suppressed by a factor q 2 /m|, where q 2 
is the momentum transfer. Because of this suppression we do not expect new contributions 
to B — > X s vv. On the other hand, the process B —>■ X s £ + £~ receives new contributions from 
the gluino-mediated 7 penguin operator and from the electromagnetic dipole operator, which 
has here a sign opposite to the Standard Model result. This change of sign implies important 
modifications both in the rate and the lepton asymmetries of B — > X s £ + £~ f28|| . The new 
contribution to the 7 penguin is however rather modest. We find that LL transitions give 

C e = ^j:U* b U ts f(m 2 /m 2 ) (39) 

9 i=l 

_ 2x(llx 2 - 7x + 2) _ 12x 4 logx 
nX) ~ 81(x-l)3 81(x-l)4' m 

where C e is the Wilson coefficient at the weak scale of the operator 

O e = ^{s a b a ) v _ A {ee) v . (41) 

Using the mass-insertion approximation and taking for simplicity m~/m 2 = 1, eq. (p9|) can be 
rewritten as 



(JSM 



5 x IQ- 2 r g . (42) 



4 Conclusions 



In conclusion we have discussed how large BR(b — * sg) can reconcile the predictions for BRsl 
and n c with present experimental measurements. Constraints from BR(B — > X s 7) require a 
precise correlation between new physics effects in the coefficients of the chromomagnetic and 
electromagnetic dipole operators. On the other hand BR(B — > X s <p) does not set any signifi- 
cant constraints on the parameter space, mainly because of the large theoretical uncertainties 
involved in the calculation. We have also shown that BR(B — > X s <fi) does not provide a good 
probe of new physics effects because the relevant coefficient C^rrib) is rather insensitive to the 
initial conditions at the weak scale of the gluon-penguin operators. 
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Finally we have discussed how a supersymmetric model with gluino-mediated flavour viola- 
tions can account for a large value of BR(b — ► sg), consistently with all other constraints from 
AB = 1 and AB = 2 FCNC processes. The theory has a rather precise correlation between 
the predictions for r g and r 7 , which is fairly independent of the specific model assumptions. Al- 
though the enhancement of BR(b — > sg) is achieved only for particular values of the parameter 
space and it is evidently not a general consequence of the model, it is interesting to know that 
supersymmetry can potentially predict BR{b — > sg) ~ 5-10%, and solve the BRsL-n c puzzle. 
Future experiments will certainly be able to test this scenario. 

Note added 

After completion of this work, we received a note by A. Kagan who informed us that their new 
results on B — > X s (j) agree with those presented in this paper . 
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Appendix 



In this appendix we report on the result of our calculation of BR{b — > s<p). As in refs. f20 , 21 
we assume a pure two-body decay and use the factorization method to estimate the relevant 
matrix elements. We have computed both the contributions of Qn and Qi2- The amplitude 
we obtain is 

A {b -> s<j>) = -^£v tb V*ig^ Uc^STLb + asi ^r bN 2N^ Cll{mb) l m ^ Lb 

( „ ™ / ml + 2Mj - m 2 A n , m b m s „ ' 



a e m b 



2 

mi 



9m b sYLb - ( 10 - ) ppRb - m s (1 
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m 2 + 2M 2 -m 2 ; 
Ml 



(43) 



where k 2 = (ml + m 2 — Mj/2)/2, C^m?,) has been introduced in eq. (|20|) and ^ is defined by 
(01^7^510) = ip^e^. This result does not fully agree with those of refs. |20], |2"Tfp|. 



From eq. (|43"D, we obtain the branching ratio 



BR(b 



T(b -> s(/>) 
F(b — ► cev) 



BRT P 



W tb V ts \ 



* 12 



SL 



\Vcb\ 2 g(m c /m b )^l(m c /m b ,m b )mf 



A(m 2 ,M 2 ,m 2 ) 



2U/2 



I^K)| 2 (2 + 2^-4-f + 2 



mf ,m 2 , m 2 N 

2^2 T ™2 ^^2 



m 2 il^ 



a s (m b )2Re(C (m 6 ))Ci 1 (m f) ) / iV 2 



57T 



m„ 



A' 2 



2iV 2 



-5mf - 20- 



2 2 

mim 



b'^s 



m 



M 2 



+ 4-JL + M1 



M; 



+ 16^-23m^| 



a e 2Re (C^,(m b )) Ci 2 (m 6 ) 
3tt A; 2 



2 2 



19mt _ + 4 



m£ 
M 2 



-23M 2 + 16^§ + m 2 ] + 



'a s (m 6 )Cii(TO fe )\ / N 2 



8ttA; 2 



2N 2 



^rn 2 M 2 



2 fi 

m£m° s 



Mi 



+ 16- 



2 4 

,m b m s 



hAm 2 m 2 — 19m. 



M 2 



4 4 

S Mi 



+ 4- 



,m£mj 
Mi 



ml 
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9 



m". 



+8^% + —M 2 m 2 - -Mi + 4^4 - 8^4 + mt | + 
M 2 2 * s 2 * Mi Ml sl 



' a e C 12 (m b ) 



87lk 2 



177 



-m 2 M 2 



2 fi 

m^m" 
Mi 



+ 112 



2 4 

. m b m s 

mT 

2 2 



- 62m?m 2 + 93m? - 4 



m 



+8 M + f M ^-f M * + 4 Mi 



M; 



a g (mb)Cii(m b ) a e Ci 2 {m b ) ( N 2 - 1 
8irk 2 8rrk 2 [ 2N 2 



4 4 

m%m s 

2 j\ /f2 



-96- 



4 2 

mim 



b" b s 



-AlmiM; - 8 



M 2 



2 fi 

^im" 
' ~Mi 



+ 4- 



,m£m, 



+ 192- 



2 4 

mim s 



M 2 



-212m 2 m 2 + 10m^ - 8 : 



4 4 

mim 



Mi 



- 128 



4 2 

mpn s 



Ml 



+ 



mYm 

) 



_- + W £|-6M<X 



+15 ^ + 8 M- 48 M + 90m -, 



(44) 



where iV is the number of colours, A(mi,m2,m 3 ) = (1 — m 2 jm\ — m 3 /mi) 2 — Am 2 m 3 /m 



1; 



2q s (m) 
3tt 



k 2 + 8z 6 - 2 8 



7T 



24z 4 ln(2:) is the phase-space correction and Vt(z,fi) 
— ) (1 — z) 2 + | is the QCD correction to the semileptonic decay rate. 



The cofficients Oj of eq. ( |i~9"| ) can be immediately obtained from eq. 



3 Numerically our branching ratio is ~ 25% smaller, including the effect of the residual scheme dependence. 
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